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We discuss the relationship between the large order behavior of the perturba- 
tive series for the average plaquette in pure gauge theory and singularities in 
the complex coupling plane. We discuss simple extrapolations of the large order 
for this series. We point out that when these extrapolated series are subtracted 
from the Monte Carlo data, one obtains (naive) estimates of the gluon conden- 
sate that are significantly larger than values commonly used in the continuum 
for phenomelogical purpose. We present numerical results concerning the zeros 
of the partition function in the complex coupling plane (Fisher's zeros). We 
report recent attempts to solve this problem using the density of states. We 
show that weak and strong coupling expansions for the density of states agree 
surprisingly well for values of the action relevant for the crossover regime. 
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1. Motivations 

Perturbation theory has played a major role in the establishment of the 
standard model of electro- weak and strong interactions. However, it is well 
known^ that perturbative series have a zero radius of convergence and that 
at some critical order, adding more terms does not improve the accuracy 
of the result. 

There exists a connection between large field configurations and the 
large order of perturbative series^ that can be illustrated with the simple 
integral 

dcjye-'^^-^^^Y.-!^ d^e-^^V^' (1) 

-OO Q • ^ — CXD 

The truncation of the exponential at order I is justified if the argument 
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is much smaller than I. However, the peak of the integrand is located at 
= a/I/. For this value of 0, the argument of the exponential is X16l^, 
which for I large enough will be larger than I. In other words, the peak of 
the integrand of the r.h.s. moves too fast when the order increases. On the 
other hand, if we introduce a field cutoff, the peak moves outside of the 
integration range and 

/+4>max , , , °° (—\V r + 't'max , , 

rf^e-i*=-A0* = ^ / d<^e-^^V'' (2) 

-4>max J —4>max 

This example suggests that one should use perturbation theory to treat 
small quantum fluctuations and semi-classical methods to treat the large 
field configurations. 

For QCD, an important challenge is to describe the large distance be- 
havior of the theory in terms of degrees of freedom which are weakly coupled 
at short distance. This question can be addressed in the framework of the 
lattice formulation. We consider the simplest case of Wilson's action which 
is the sum over the plaquettes in the fundamental SU{N) representation: 

5 = ^(1 - il/N)ReTriUp)) . (3) 
p 

With the usual notation (3 = 2N/g'^, the partition function reads 

Z = l[ [ . (4) 

For iV = 3, this theory has no phase transition when one goes from small 
coupling (large /?) to large coupling (small /3). Recently, convincing argu- 
ment have been given^'^ in favor of the smoothness of the RG flows between 
the two corresponding flxed points. Consequently, there does not seem to 
exist any fundamental obstruction to match the two regimes. One general 
question that we would like to address is if it is possible to construct a 
modified weak coupling expansion that could bridge the gap to the strong 
coupling regime. 

Lattice gauge theories with a compact group have a build-in large field 
cutoff: the group elements associated with the links are integrated with dUi, 
the compact Haar measure. For SU{2) and SU{3), the action density has 
an upper bound which is saturated when the group element is a nontrivial 
element of the center. It is remarkable that this formulation has a UV reg- 
ularization and a large field regularization that preserves gauge invariance. 
Does the presence of large field cutoff means that pcrturbative series are 
convergent? Not necessarilly, because in constructing perturbative series. 
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one dccompactifies^ the gauge field integration which at low order amounts 
to neglect exponentially small tails of integration, but modifies the large 
order drastically. 

In continuous field theories, it is expected^"^ that the large order of per- 
turbative series can be calculated from classical solutions at small negative 
A for scalar theories or small negative g"^ for gauge theories. For lattice gauge 
theory with compact groups, the theory is mathematically well defined at 
negative ,i. e., negative /3, but some Wilson loops oscillate when their 
size change and the average plaquette has a discontinuity when g"^ changes 
sign.^° For small negative g"^ , i. e., large negative /3, the behavior of expec- 
tation values is dominated by the large field configurations. This statement 
can be made more precise by introducing the spectral decomposition 



with n{S) the so-called density of states. It is clear that for large negative 
/?, what matters is the behavior of n{S) near Smax- 

We can make the discussion more concrete by considering the case of the 
single plaquette model-^^ with SU{2) gauge group. In this case, Smax = 2 
and we have 



The weak coupling expansion of Z is determined by the behavior of near 
S = 2. When we expand V2 — S about S* = 0, the large order of this 
series is determined by the cut at S* = 2. After integration (from to +00) 
over S, the series with a finite radius of convergence becomes asymptotic 
(inverse Borel transform). Alternatively, we can maintain the finite range of 
integration and construct a converging weak coupling expansion but with 
coupling dependent coefficients as done by Hadamard a century ago in his 
study of Bessel functions. It would be interesting to know if the features of 
the one plaquette model persist in the infinite volume limit of lattice gauge 
theory. 

In these proceedings, we discuss the large order behavior of the weak 
coupling expansion of the plaquette (Sec. 2) and the possibility of defining 
its non-perturbative part (the glnon condensate? Sec. 3) . More details can 
be found in Ref.^^ In Sec. 4, we present numerical results concerning the 
zeros of the partition function in the complex coupling plane^^'^^ (Fisher's 
zeros). This problem could be solved using the density of states. Prelimi- 
nary numerical results concerning the density of states are provided in Sec. 





n^pi.{S) = -^/S{2-S) . 



(6) 
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5 where we also show that weak and strong couphng expansions for the 
density of states agree surprisingly well for values of the actions relevant 
for the crossover regime. After the conference, we wrote a more detailed 
preprint^^ concerning this question. 

2. Perturbative series in lattice gauge theory 

In this section, we denote Np = L^D{D — l)/2 the number of plaquettes 
and the average plaquette: 

P(/3) = - {l/N)ReTr{Up))j . (7) 

The weak coupling expansion of this quantity 

P{P):^J2^rar"' + --.. (8) 
m=l 

has been calculated up to order tcn-'^^ and 16."'^'^ Series analysis"'^*'^^ suggest 
a singularity P oc (1/5.74 — l//3)^ °* and consequently a finite radius of 
convergence. This is not expected since the plaquette changes discontinu- 
ously^° at /? — > ±oo. This behavior is also not seen in the 2d derivative of 
P and would require masslcss glucballs. 

A simple alternative is that the critical point in the fundamental-adjoint 
plane^" has mean field exponents and in particular a = Q. On the Padj. = 
line, we assume an approximate logarithmic behavior (mean field) 

-dP/dp oc H{i/pm - i/Pf + r') . (9) 

r cannot be too small (absence of singularities) or too large (this would 

create visible oscUations in the perturbative coefficients). From these con- 
straints, we get the approximate bounds^^ 0.001 < F < 0.01. 
Integrating, we get the approximate form 

hr" ^ C(Li2(/3- V(/3~i + iT j) + h.c , (10) 

k=0 

with 

Li2(ar) = 5^a;Vfc' • (11) 

fe=0 

We fixed F = 0.003 and obtained C = 0.0654 and /3^=5.787 using of 69 
and 610. The low order coefficients depend very little on F (when F < 

0.01), larger scries arc needed to get a better estimate of F for SU{3). 
It interesting to notice that we get very good predictions of the values 
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Table 1. 


Predicted values of bm 


with the dilog 


model. 








order 


predicted 


numerical'^® rel. error 


1 


0.7567 


2 


-0.62 


2 


1.094 


1.2208 


-0.10 


3 


2.811 


2.961 


-0.05 


4 


9.138 


9.417 


-0.03 


5 


33.79 


34.39 


-0.017 


6 


135.5 


136.8 


-0.009 


7 


575.1 


577.4 


-0.004 


8 


2541 


2545 


-0.0016 


9 


exact 


11590 




10 


exact 


54160 





of bs,b7,. . .1 Wc believe that these regularities should have a Feynman 
diagram interpretation. 

Another possible modeP^'^^ is based on IR renormalons 

J2 bkr" ^ K r dte-^' (1 - t 33/167r2)-i-204/i2i (12) 

/3 = /3(l + rfi//3+...) (13) 

ti — corresponds to the UV cutoff while t2 ~ 167r^/33 corresponds to 
the Landau pole. For <2 = 00 we get a perturbative series with a factorial 
growth, in contrast with the previous model which had a power growth. Un- 
fortunately a clear distinction between the two types of large order behavior 
requires numerical calculations at order larger than 20. 



3. The Gluon Condensate 

Using the two large order extrapolations described in the previous section, 
one can see good evidence^^ for 

P{P) - Ppert.iP) ^ B{a/ro)\ (14) 

with a(/3) defined with the so-called force scale,^^'^^ and Ppert appropriately 
truncated for the second model. B is sensitive to resummation. B ~ 0.7 with 
the bare series^^ and 0.4 with the tadpole improved series.^'' 

It is tempting but potentially dangerous to try to relate B to the nu- 
merical value of the gluon condensate. If we identify^^ for AT = 3 

P(/3) - Ppert. m < -GG >, (15) 
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wc obtain for tq = 0.5 fm and B = 0.4 that < f GG >~ 0.035 GeV"^ which 
is about 3 times the original estimate. It is not clear to us that there is a 
precise correspondence between the continuum and the lattice definitions. 
Also, different values for the continuum value have been proposed. A nega- 
tive value, subsequently criticized^^ can even be found. For these reasons, 
some authors^^ prefer to use a value with error bars when estimating ag 
at some large scale. If we use the correspondence between the lattice and 
the continuum discussed above at face value with B = 0.4 — 0.7, these error 
bars should be multiplied by a factor 3 to 6. 

4. Zeros of the Partition Function 

The existence of complex singularities near the real axis can be tested by 
studying the complex zeros of the partition function.^" The basic technique 
is the reweighting^^'^^ of action distributions at given f3o 

Z(/3o + A/3) = ^(/3o) < exp(- A/35) >/3„ . (16) 

As mentioned in the introduction, Z{6) is the Laplace transform of the 
density of states n{S). For SU{2) with even nunil)crs of sites in each direc- 
tionio Z(-/3) = e^l^^^Z{l3) and n{S) = n{'2M'p - S). In the crossover, we 
have attempted local parametrizations:^^'^^ 

n{S) a Q-{aiS+a,S-+a,S-+a,S-) (^7) 

For SU (3) on a 4^ lattice, stronger deviations from the Gaussian behavior 

are observed than for SU{2). This is illustrated in Fig. 4 that has a larger 
scale than its SU (2) counterpart. The histogram were made with 50,000 
configurations prepared for a study of the third and fourth moments. 
As the volume increases, these features tend to disappear in the statistical 
noise. These local model give results that can be compared to MC rewcight- 
ing in the region where the errors on the zero level curves for the real and 
imaginary parts are not too large. An important consistency test is to 
show that approximately the same zeros are obtained for different (3q . This 
is illustrated in Fig. 2. We plan to pursue this study using the more global 
information provided by the density of states. 
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SU(3); 6x6x6x6; beta =5.58 



binwidth: 7.0 
50000 contigs. 




3400 3500 3600 3700 3800 3900 4000 4100 

S 



Fig. 1. The relative difference with a Gaussian distribution for 50,000 values of 5 in an 
histogram with 100 bins for a SU(3) pure gauge theory on 4** and 6* lattices at /? = 5.58 
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SU(3); 4x4x4x4; /? =5.54 d=0.25 

0.20 



0.15 



0.05 



0.00 

5.40 5.46 5.51 5.57 5.62 5.68 

Re. 13 



SU(3); 4x4x4x4; /? =5.55 d=0.25 

0.20 



0.15 



0.05 



0.00 

5.41 5.47 5.52 5.58 5.63 5.69 

Re. (3 



Fig. 2. Zeros of the real (circles) and imaginary (crosses) for 5(7(3) on a 4'* lattice 
using reweighting of MC data at /3o = 5.54 and 5.55. The small dots are the values for 
the real (green) and imaginary (blue) parts obtained from the 4 parameter model. The 
crossed box at ((5.541,0.103) above and (5.542,0.102) below) have been obtained with a 
perturbative method. Red boxes: boundary of the MC confidence region. The solid line 
is the circle of confidence of the Gaussian approximation. The locations of the complex 
zeros are consistent. 
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5. Approximate form of the density of state n{S) for 
SU{2) 

We assume the following form: 

n{S) = e^^f^'^/^^^ . (18) 

In the infinite volume limit, f{x) becomes volume independent and can 
be interpreted as a (color) entropy density. In the same limit, we have the 
saddle point equation 

fix) = P . (19) 

This is the analog of the familiar thermodynamical relation dS/dE = 1 /T. 
Knowing f{x) amounts to solve the theory (in a thermodynamical sense). 
For SU{2), f{x) = /(2 — x) (symmetric about 1) and we don't need to cal- 
culate / for a; > 1. We have constructed weak coupling and strong coupling 
expansions for SU (2) and compared with numerical data on a 4^ lattice. 

For the weak coupling, wc use the large beta expansion near x = 0, 
X =< S/Afp >= 1^ + 0.156^ + . . . . We assume f{x) ~ A ln{x) + B + 
Cx + . . . , plug the expansion in the saddle point equation and solve for A 
and other unknowns. We obtained f{x) = jln{x) + 0.208a; + 0.0804x2 + 
For the strong coupling, we use the low beta expansion^^ to solve 

4 96 



near a; = 1 for SU{2), and obtained a; - 1 =< S/Np - 1 >- ^ 



lik + Mmo + ■■■ + 68i9i30659Sioo ' ^ith periodic boundary conditions, the 
coefficients have no volume dependence for graphs with trivial topology 
(volume dependence should appear at order /J^^). Solving the saddle point 
equation for an expansion about 1, we get f{x) = —2{x — 1)^ — |(a; — 1)** + 

I 163150033(z-l)^° , 

255150 -|- . . . . 

The numerical construction of n{S) by patching was done by A. Den- 
blcykcr and is illustrated in Fig. 3. Scries expansions of f{x) are compared 
with the numerical data in Fig. 4. Note the good agreement in the inter- 
mediate region. After the conference, we wrote a more detailed preprint^^ 
where details can be found. 
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0.0275 - 
0.0300 - 
0.0325 - 
0.0350 ■ 
0.0400 
0.0450 
0.0500 - 
0.0600 
0.0700 
,0.0800 - 



1 400 1 420 1 440 1 460 1 480 1 500 1 520 1 540 1 560 1 580 1 600 1 620 



200 400 600 



1000 1200 1400 1600 18( 



Fig. 3. n{S) oc Pi3^^{S)e^^^ for different /3o (patching). Collection of overlapping data 
(A. Denbleyker). 



6. Conclusions 

The density of state show a nice overlapping of the strong and weak coupHng 
expansions. We plan to use the density of states to study the Fisher's zeros. 
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SU(2) 4x4x4x4 




500 1000 1500 

S 



Fig. 4. ln{n{S)) numerical (black), strong coupling at order 16 (blue) and weak coupling 
with dilog model at order 10 (green) for SU (2) on a 4'* lattice. 

We need numerical confirmation of guesses made for the weak coupHng 
expansion for SU{2) where Im(3F / RejSp is 5 times larger than for SU{'i) . 
We need better understanding of the lattice and the continuum definitions 
of the gluon condensate. We need a better understanding of the large order 
behavior of QCD series in terms of the behavior at small complex coupling 
(a picture analog to metastability at A < for the anharmonic oscillator^"^). 

This research was supported in part by the Department of Energy under 
Contract No. FG02-91ER40664. 
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